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I , Abstract 

}—( ■ 
-i—i , We propose a new approach to the self-consistency equation, which arises in the problem of the 

- motion of a hole in a quantum antiferromagnet, appropriate to the case of small exchange energy 

"_{ " J. The functional equation for the Green function is transformed into a differential equation; its 

solutions are analyzed and compared to the existing numerical calculations. This method allows 

one to study the limit of J — > 0. Application to other strongly correlated electron systems is 

1-^ discussed. 

a 

: 1 Introduction 

The interplay between doping and antiferromagnetism is believed to play an important role in high Tc 

^ , superconductors. Hence the motion of holes in an antiferromagnetic background has received much 

f~^ ' attention in the last two decades. After the pioneering work of Brinkman and Rice [I], where the 

0^ , Hubbard model with t/ = oo is considered, many publications on this subject are known [5], [3, 01, 

|S]. Often this problem is faced by applying the model t-J Hamiltonian, acting on the space with no 

doubly occupied sites, which is obtained from the Hubbard large U model. One may use, for instance, 
\^ • the Holstein-Primakoff transformation (see jj) to rewrite the Hamiltonian it terms of interacting 

f^ , spinless holes and spin waves in order to make it similar to that of the polaron problem. Then in the 

self-consistent Born approximation the selfenergy may be calculated ([2], [3) 



s(k, Lu)^Y. /C^' q)G('' - q> ^ - E,) 

G , and the result is an integral equation for the Green's function of the form 

^' G{k,u;)^ ^ », .i,., TTT (1) 

r> ' where /(fc, q) contains information about the coupling of the holes to the spin excitations and is of the 

j^ ■ order of t , t being the hopping rate of the hole; Eg is the energy of the spin excitations, which is of 

order of J, the antiferromagnetic exchange energy. This self-consistent integral equation contains all 
the information about the hole spectrum ; numerical analysis has shown, that it consists of a broad 
band of incoherent states and near to its bottom - a peak in the density of states (DOS), corresponding 
to a quasiparticle state (see e.g.|21)- 

It is often claimed that the k-dependence is not too important , [3 ■ Then it is worth to analyze 
thoroughly the k-independent case ( the Ising limit), when |^ take a simpler form: 

G(C^) = T^, y-y (2) 
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(We have absorbed a factor z (the coordination number) into i^) 

This functional equation (jSJ is investigated in the Uterature in relation with the polaron problem 
|3|, |H], the problem of the hole motion, described above, but also in other physical problems. We 
mention another example here - the one-dimensional anisotropic ferromagnetic chain, which is studied 
in relation with experiments on the magnetic salt CoCli-'2.il20 |S], |1U| . 

The resulting DOS consists of a series of sharp delta functions, as shown i.e in (S| by numerical 
calculations. This picture is in contrast with the smooth DOS of an incoherent band for J = found 
in and it is not clear how the limit J ^ can be obtained. 

In order to understand the behavior of the system it is useful to realize, that the motion of the hole 
creates "strings" of overturned spins; their field confine the hole, thus leading to selftrapped states 
centered at the original hole position ^21) d]i dji |15| . 

Here we present an approximated treatment of eq. (O suitable for small values of J and which 
allows us to consider the limit J ^ 0. We transform the functional equation into a differential 
equation and we show that the solutions of the latter present a picture of the DOS in agreement with 
the existing numerical calculations. 

The paper is organized as follows. First we discuss the existing solutions of eq. |5J, taking 
as a model the investigations on magnetic salts. Then we transform the functional equation into a 
differential equation and we solve it. At the end we analyze the resulting DOS and outline the changes 
induced by the interaction J. We show that our solutions are in line with the previous calculations. 

2 Discussion of the existing analytical solutions 

In order to outline the simplicity of the method, proposed here, we first recall the exact solution, 
following in the main the calculation for the case of magnetic salts. Pj, |11| . 

First of all notice that the functional equation Q is equivalent to the continued fraction : 

1 
G= 






UJ — J — 



UJ-2J - ... 
which in turn can be thought of as solution of the series of difference equations: 

c„(w - nJ) + t(c„+i + c„_i) = (3) 

for n = 0, 1 , 2 , . . . Indeed, from ^ one obtains 



Co 



t^ 



oj — J - 



UJ-2J - ... 

The spectrum may be found from the boundary condition c_i =0. 

Compare now to the case of magnetic salts |11| . The Heisenberg Hamiltonian with a magnetic 
field Hq is simplified to reflect the main properties of the system. So, one starts from the eigenstates 
of the Ising: 

N N 

H' = -2r- ^ StSU, + iHo Y, St , 

i=l i=l 

which consist of single clusters of n-adjacent spin deviations with respect to the ferromagnetic ground 
state, with selfenergies 

4")=2J-+n7iJo, (4) 

and constructs Bloch functions. In this basis is calculated the perturbation due to the anisotropic 
term, which is the most important in the case of ferromagnetic salts: 

N 



jja ^_ja ^(5+5-+ ^ + H.C.) = HP'^'-K 



For the coefficients in the perturbation series one obtains the same type of coupled equations as Q the 
only difference being due to the fact, that the perturbation TJp^''* couples states with spin deviation 
differing by 2 0: 

wyHo - J° cos(fca)(c„+2 + Cn-2) = (E ~ 2J^)c„ 
with a - the interatomic distance and — — < fc < — . 

a — — a 

The correspondence between the two cases, discussed above, can be expressed as follows: 

(E ~ 2 J^) -^ uj; jHo -^ J; J" cos(fca) -^ t. 

For the simplest case J = (or in the absence of magnetic field, when one considers the magnetic 
salts) (jSl may be solved taking — — = — ^ — T. The "transfer matrix" T obeys the equation: 

Cn ^n — 1 

Tuj + t{T'^ + 1) = 
and we have 



T 



(-U + ^/lo^ - AtA /2t 



We obtain continuous spectrum in the region ~2t < lo < 2t . This corresponds to the semispherical 
DOS found in p. Let us put w = 2tcosA, where A e (0, ^). For T we have two solutions: T = e^*"^, 
from which we can construct the expressions for c„, corresponding to the appropriate boundary 
condition : c„ = ^sin((n + 1)A). 

Now, let us consider the solution of JSJ in the general case. When we realize that the same 
recurrence relation satisfy the Bessel's functions, we can write this solution at once: 

c„ = ^J(._„)/2(^). (5) 



The boundary condition is then 



'^(t+o/Ht'-'^- 



The spectrum can be found as the zeros of the Bessel functions are tabulated (see i.e. jl6|). 

It is possible to obtain more transparent expressions using the appropriate double asymptotic 
expansion of the Bessel function in both order and argument [Hj , |17|: 

Ja{h^Ji) « \ cos[^(tana - a) - 7r/4] , +0(^"^/^) 

V TT/i tan a 

where 6=1/ cos a, < a < 7r/2. 

Setting LO — 2t cos A as earlier, we have 



2t 

LO 



COS A cos a 

i.e A = a and the zeros of the Bessel functions can be calculated from the condition: 

3 

sin A — A cos A = ti{v H — )J/2t , 

where i^ = 0, 1, 2.... From this formula the spectrum can be calculated numerically. 

Another type of approximation to the recurrence equation© consists of considering the coefficients 
c„ as slowly varying functions of n, which is justified in the limit of large n and is on the same lines 
as the calculation in ^^l- Replacing c„ by a continuous function c(n) and expanding c{n± 1) about 
n, we obtain the following equation: 

c{n){io ~ nJ - 2t) + tc" {n) ^ Q (6) 



which may be regarded as a one dimensional Schrodinger equation, describing the motion of a particle 
of mass l/2i in the linear ("wedge") potential U{n) — nJ. The boundary condition c(— 1) = 
is equivalent to the introduction of an infinite potential wall. This problem is well known to yield 
solutions in terms of Airy- functions. In WKB approximation the Bohr-Sommerfeld condition reads: 



p{x)dx^ ^/{uj - 2i + J - xJ)/t dx = 7r(j/ + 3/4), i^ = 0,l,2.... 



where a; = n + 1 , x^ = [lo — 2t + J) / J and one obtains the spectrum 

io-2t + J^ (37r(i. + 3/4)/2)2/3ii/3 j2/3^ 

The spectrum, found here, consisting of discrete levels, should merge to form a continuous band, 
when J — > 0. But this fact cannot be expressed analytically, the limit J — > being singular. 

3 A new approach to the functional equation 

Let us write the equation for the self-energy, corresponding to eq. |(2J): 



For small J we can expand the function S : 





^\ 


^1 T 

UJ — J 


- [E(cj) - JS' 


H + 


i(J)2l]' 


'H + -] 


or 




I](w- J- 


S) + JS'S - 


lUf 


■Y."T. + . 


..^e 


When J — Q the resulting equation 
has the well known solution 


S(tj-S) 


= f 







(8) 
(9) 



{uj ± ^^^2 _ 4i2)/2 (10) 

which gives a semispherical density of states (DOS) in the region \uo\ < 2t. 

The differential equation cannot be integrated analytically. Nor is it possible to treat the terms 
~ J as a perturbation; so, in order to find an appropriate solution we proceed as follows: 

We will transform first into a higher order equation, 

I]'(co - J - S) + E(l - S') + J(S')' + JS"S - i(J)2E"l]' - .... = , (11) 

having in mind, that only part of the solutions of the latter are also solutions of ©. Then we have 
to choose the integration constants as to reproduce the correct result for the equation of lower order. 
We are interested in the solution for the energy region of the quasi-band states, as we want to 
know how the DOS, generated from H10() . changes, when one includes a small interaction J. We look 
for a solution, which is sufficiently smooth ^ . So we neglect the terms ~ JS" and obtain the following 
equation (x = lu — J and y = S): 

Jy'^ + y^i2y-x)y'. (12) 

We believe this equation contains the essential part of the information needed to analyse the change 
of the DOS in the band. In the Appendix we outline the main steps of the calculation. The solution 
is obtained in parametric form: 



J^-A/^^ + i^-i/^^-t-/^ (13) 



^Because of this approximation details of the DOS structure near the band edges are lost 



2P-1 l^_P 



P 



p-1 



J J 2p~l 
2"" 2" p 



l~p 



■ arctan 



l-p' 



where D is an integration constant. In order to reproduce the known resuhs without J, we have to 
choose D — t^. Indeed, let us take J = 0, we have then: 



2 n P 
y ^ D- 



ana x = y 



p~ 1 p 

whereby, getting off the parameter p , we obtain the direct relation between x and y: 



D 

X = \-y 

y 



and we have again pU|) . 

Another useful relation reads 



P=^[22/-2;t V(2y-a;) 



Uy]. 



(14) 



(The sign is to be chosen properly, as to reproduce at J ^ the starting equation y' — y/{2y — x) ). 

The next step is to transform these relations in some formulae, which enables us to analyze the 
change in the DOS under the influence of J. 

We rewrite the second equation in H13|l in the form: 



J 2p~l 



p-1 



where 



J I V 

t = y D — i— arctan . / 

2 V l-p 



(remember, that vD = t). Then we introduce a parameter q: 



(2p - 1)2 _ £2 ^ 
p{p - 1) " F " "'^ 



and find, making use of (|14|l : 



y 



. + J+jMi±il±^ + 2J. + J2 



9(4 + q) 
Then the parameter i is obtained from the equation 



i[ 1 + 



t = t+ -Arth ^^-= 

2 2i 



and the solution for y is: 



y 



+ y/i^ - 4i2 



J 



A + q 



2 y/i^-Ai\ 



4 + g 



(15) 



(16) 



It is clear, that the band width is now determined by the value of t, which is function of x to be 
calculated from H15|) . This transcendental equation can be solved numerically, but several properties 
may be revealed by analytical study of the results. 




Figure 1: Plot of the function f{i) (the symbol t on the plot stays for i and the hopping rate t = I) 
for 3 different values of J :0,2 ;0,4 and 0,6. 

4 Analysis of the solution 

We will still use for convenience the variable x — oj — J/2 (a shifted energy). From Hlt)|) we have the 
self-energy 



i+^/S;^- 4P / J 1 

^ — T^ I 1 ~l~ 77" 



2 ^/S:^ - 4? 



(17) 



and we can calculate the Green's function 



G = 



i + j-Y. 



in terms of the parameter t: 



G = 



(^2 - AP){i + ^) + v/i2-4p(i2 _ 4P + ^) 
2i2(:j2_4p_^) 



(18) 



The parameter t plays in this formula a role similar to that of the hopping rate t in the absence of the 
interaction J. In the latter case t determines the width of the incoherent band (4t). Now t depends 
on X (on the energy) and it is not obvious how the band width changes. 
The relation to determine t H15|l may be rewritten in the form 



i = 2i coth ( -{i~t) ] = /(t), 



(19) 



so, the function t(x) can be found graphically using the plot of /(f) (Figure^. 



First we consider the region above the top of the band. For high energies (large i > 0) sohition 
is obtained from the curve f{t) above the jump at i = t{t = 1) and we have i > 1. The physically 
meaningful solution for t{x) is obtained on the left branch of the curve , where t{x) ^ t a.s x ^ oo. 

Once t is known, one can calculate the Green's function. In this region of x our solution satisfies 
X > 2t and G is real. We may use the following expression for it: 



G 



tank {^{i-t)) 



i+i 



sinh (j(t-t)) 



(20) 



The function/ (i) - (|19ll - has a minimum in this region and the corresponding value of t, t„i is to 
be evaluated from the relation : 

sinh((8/J)(i™-i)) =(8/J)t™. 

Below the minimum of / real solutions of the transcendental equation 1)19(1 do not exist, so it 
represents the upper limit (the top) of the band. We have calculated the minimum value of / (or the 

top of the band Xtop) for 3 different values of the exchange interaction J( see Table 1). 

Table 1 



J/t 


0,2 


0,4 


0,8 


Xtop/t 


2,275 


2,489 


2,8637 


ttop/t 


1,11222 


1,1934 


1,3285 


ig/t 


-1,649 


-1,40 


-1,23 



We notice the almost linear increase of the top of the band with J beginning from its value of 2t 
at J = 0. 

For X < fmin wc havc to calculate the complex t = ti + it2 and then find the complex Green's 
function which gives the spectral function (and the DOS in the band). From ((191 ) we obtain the 
following system for the real and imaginary part of t: 



— X 

J 



277 



S; 



sin(2?7) 
2,7 _2^+f 



sin(277) sin(2^) 



(21) 
(22) 



where the following notations are used: 

A 

B = cosh(20 +cos(277). 

In terms of these variables we may write the Green's function in the form: 

cosh(^ — iiTi) 



4 4 4 



^" 



[cosh(^ — ii]) — B] 



1- 



2?7 



sin{2ri) 



cosh(^ — irj) 



(23) 



We study analytically G in the center of the unperturbed band (x — 0) and in its vicinity. For 
a; = 0, ^ = and 2rj = n . Let us define 

2^ ^5; 2rj^TT + A. 

and consider the case d ^ 1 and A <C 1. From ((21(1 and 1(22(1 we find 

X . Sit 



1 + [Jn/sy 



26 + 8/ J 



(24) 



and for the Green's function 

2 — TT — ITT—- 
A 

As (5/A « —8/Jn roughly, we obtain for the real and imaginary part of G (the latter gives the DOS) 
in the vicinity of a: = 

ReG^ W8)(2--)^, j,,G^ ^ ,. (26) 

1 + [(2 - tt) J/8]' 1 + [(2 - tt) J/8]' 

Another important feature of the calculated Green's function is it's behavior in the region of the 
lower part of the band. Consider now the negative values of the energy (of x ) . At i; — > — c» we find 
again i ^ t , but i < t (on the left side of the jump of the function /(t) ), as it is expected. The ratio 
\x/2i\ > 1 , so we have real values of the Green's function and we are outside the band. 

On the other hand, as sinh {j{i — t)) < for x < 0, we see from (pn|) that the Green"s function 
has a pole. Hence, we obtain a quasi-particle state. The corresponding energy is calculated from the 
condition: 

(4/ J)t = - sinh ((4/ J) (i - t)) (27) 

and making use of l|19() : 



X 



y^l + sinh'((4/J)(i-i)) 



2t sinh((4/J)(t-i)) 

We find 

,Ji^-4P{i) = J/2 (28) 

for the energy of the quasiparticle level. The calculated energies Xq of this level for 3 different values 
of the parameter J are given in Table 1. 

For negative energy smaller in absolute value we arrive at the point, where the condition 2i « x 
is satisfied. (For small J this point is just next to the quasiparticle level). The latter condition is the 
bound of the real values of coth ( j(t — t)). So , we are at the lower limit (the bottom) of the band. 

5 Conclusions 

. We developed here a different method for solving the functional equation (0) by transforming it into 
a differential equation. This approach may be useful when the exchange interaction J is not large. 
Under this approximation the region of sharp peaks in the DOS merge into a smooth band. We can 
see how the top and the bottom of the band change with J - they are shifted to higher energies, 
numerical results obtained for the more general equation iQ . Similar behavior is obtained as a result 
of numerical analysis of the more general equation |^(see e.g. |H|,|3]). 

Also a bound level (quasiparticle state) immediately under the incoherent band is found, in agree- 
ment with previous numerical calculations j^. 

So we have a picture of the DOS, which is very close to that known from the literature. 

We also obtained an explicit expression for the Green's function near the center of the band. 

We believe, that the proposed method of treating the functional equation can be used as a com- 
plementary tool to study the DOS, generated by the self-consistency equation jSJ. 

Numerical calculations of the DOS in the entire band are in course. 

6 Appendix 

Here we write down the main steps in finding the solution of H12|) of the text. 
First we label y' = p and we write it in the form: 

x^2y-Jp~-. (29) 

P 



We then differentiate to obtain: 



Tfien we are at tfie equation 
and its solution 

y{p) = 

Returning to (|29|l we have: 

x{p) ■ 



dx = 2dy — Jdp dy + -hrdp = — 

p p'' p 



dy y ' J P 



dp 2pil-p) 2(1 -p) 

(arctanJ— Vp(f-p)j 



l-p 



--I 



2p-i r~p 



i-p 



C-i.rct..~ 



2 \ l~P 



J 
2' 



Redefining the integration constant, we can write C ,/ ^r- — \h^~^ ^^d then system 



y^p) - a/^^ + {p-^2 Vt5"'^'*""v^ ^^'^^ 



2p-i /n p J J 2p~i r-f- ry 

X = . / U : — — — — 1 / arctan . 



p-l 2 2 p y l-p \ l-p 

define our parametric solution of p2|) . 

Without J the solution is readily obtained, writing Hf2() as 2yy' = xy' + y and integrating: 

y^ + Ci = xy; 

here the integration constant Ci is identified as Ci = D. 
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